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Abstract 

We address the problem in Navier-Stokes isotropic turbulence of why the vorticity accu- 
mulates on thin sets such as quasi-one-dimensional tubes and quasi-two-dimensional sheets. 
Taking our motivation from the work of Ashurst, Kerstein, Kerr and Gibson, who observed 
that the vorticity vector u aligns with the intermediate eigenvector of the strain matrix S, 
we study this problem in the context of both the three-dimensional Euler and Navier-Stokes 
equations using the variables a — £• S$, and x = £ x where £ = uj/uj. This introduces the 
dynamic angle (/>(x, t) = arctan which lies between u) and Su. For the Euler equations 
a closed set of differential equations for a and x is derived in terms of the Hessian matrix 
of the pressure P — {p,y}. For the Navier-Stokes equations, the Burgers vortex and shear 
layer solutions turn out to be the Lagrangian fixed point solutions of the equivalent (a, x) 
equations with a corresponding angle <f> — 0. Under certain assumptions for more general 
flows it is shown that there is an attracting fixed point of the (a,x) equations which corre- 
sponds to positive vortex stretching and for which the cosine of the corresponding angle is 
close to unity. This indicates that near alignment is an attracting state of the system and is 
consistent with the formation of Burgers-like structures. 
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1 Introduction 



As early as 1938 G. I. Taylor |jj showed that in isotropic Navier-Stokes turbulence vortex 
stretching has a major effect on vorticity production and dissipation. One of the many inter- 
esting features of high Reynolds number turbulent flows, illustrated beautifully by modern flow 
visualization methods, is the fact that vorticity is not evenly distributed throughout a flow but 
has a tendency to accumulate on 'thin sets'. The morphology of these sets is characterized by 
the predominance of quasi one-dimensional tubes or filaments and quasi two-dimensional sheets 
although, being fractal in nature |2|, |3], |j, many individual vortical structures may be neither 
precisely one nor the other. This morphology makes it clear that the direction of vorticity as 
well as its production plays a significant role in the self-organization processes through which 
apparently intense regions of vorticity appear to metamorphose into approximate thin coherent 
structures. The challenge facing the theorist working with the three-dimensional Navier-Stokes 
equations is to explain these geometric accumulation processes and their consequences even 
though a detailed description of their precise topology is beyond our reach at present. The 
subtle fine spatial structure of these sets indicates that an alternative approach is needed to the 
conventional one of Navier-Stokes analysis where L 2 -spatial averages are taken over the whole 
domain when estimating the effect of the vortex stretching term. 

In 1951, Townsend |^, |6| indicated that something like the thin sets described above might be 
relevant in turbulent flows and in the last two decades their occurrence has been a regular theme 
of the literature @, |, g, [fij 0, (j^, ||, 0, P, £§, 0, |l|, g^, ||, |l|, ||, ||, g 



26, 27, 



They would also seem to be the visual manifestation of the vorticity alignment phenomenon first 



reported by Ashurst et al [14] who observed that in driven simulated isotropic Navier-Stokes 
turbulence with a Taylor microscale Reynolds number of 83 the vorticity vector u> has a tendency 
to align with the intermediate eigenvector of the rate of strain matrix 



S, 



ij 



2 Kj + u h i 



(1) 



This observation was based on a study of the probablility density function (PDF) of the cosine of 
the angle between their directions. This important alignment process has also been confirmed in 



other numerical simulations |T5|, 16, 18, [2(J, 26, |28[] as well as in turbulent grid flow experiments 
conducted by Tsinober et al [29, |3^, [n], 32] and it has been suggested that it is a kinematic 
process [33, |34) . In fact preferential vorticity alignment appears to be a universal feature of 
smaller scale structures in such flows even though larger scale features may vary from flow to 
flow. A common feature seems to be that vorticity is concentrated in tubes of width intermediate 
between the Taylor and Kolmogorov microscales pi, 111, pi |fl, [28| with viscous dissipation 
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occurring in annular regions around the tubes which themselves form from the roll-up of vortex 



sheets f2q| . 

In isolation from each other and in an ideal sense, vortex tubes and sheets can be thought 
of as Burgers vortices and Burgers (stretched) shear layers respectively for which there are 



known exact solutions of the Navier-Stokes equations in certain special cases |25|, 13. 24, [35|| . 
Of course the Burgers solutions are idealized because their strain fields are not linked back 
to the local vorticity field as would be the case in a real flow. Nevertheless they have the 
great merit of providing us with exact solutions which allow an interpretation of the processes 
involved. A Burgers vortex tube ideally has one positive and two negative eigenvalues of the 
strain matrix (axial strain) while a Burgers shear layer solution has two positive and one negative 
eigenvalue (biaxial strain). Ashurst et al H] reported that their eigenvalues occurred in the ratio 
3:1: —4 but, as Moffatt, Kida and Ohkitani have pointed out [25], vortex tubes can still survive 
in regions of biaxial strain provided they are strong enough. As best illustrated by the graphics 
of Vincent and Meneguzzi |26[| , vortex sheets deform neighbouring sheets, curling up like potato 
crisps (chips) most strongly when rolling up into tubes. It is significant that they found that the 
tendency towards alignment between u and the intermediate eigenvector of S occurred before 
the roll-up of sheets into tubes. Burgers vortices have been produced in the laboratory by 
Andreotti et al [^] from two experimental set-ups which allow them to study the phenomenon 
of vortex stretching in detail. While the Burgers solutions are extremely valuable to the theorist, 
it is necessary to formulate a theory which demonstrates why a flow should evolve to such states 
in the first place. 

For the three-dimensional Euler equations the situation is somewhat different and com- 



plicated by the suspicion of a finite time singularity in the vorticity field (see Kerr [37, |38||). 
Constantin, Fefferman and Majda ^0|, however, have recently shown that singularities in 
solutions of the three-dimensional incompressible Euler equations can be ruled out if the velocity 
is finite and the direction of vorticity is smooth but they cannot be ruled out if the direction of 
vorticity is not smooth. For instance, blow-up cannot be ruled out if two intense vortex tubes 
collide at a nontrivial angle. In fact, Kerr's computations predicting such a singularity used 
anti-parallel vortex tubes as initial data p7], (see also EL [42|]). The result of Constantin, 



Fefferman and Majda [39| is related to that of Beale, Kato and Majda [43| which has provided 
the main criterion for understanding the growth of vorticity in Euler flows for more than a 
decade. 

We begin by setting up the notation and then proceed to give some preliminary definitions. 
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Let us consider the incompressible three-dimensional Euler equations in vorticity form 

Du 



Dt 

and the Navier-Stokes equations 



(2) 



Du 

— = a + i/Aw, (3) 

with a divergence free velocity field divu = 0. The material derivative is defined in the usual 
manner as 

The vortex stretching vector a appearing in @ and (|3|) is given by 

o-j = ujjUi.j = SijUJj. (5) 

(j is written as <x = Su: on the right hand side of (|5|) because u> sees only the symmetric part of 
Uij. The approach taken in this paper is to consider the dynamics of the local angle </>(x, t) which 
lies between the vorticity vector u and Su at the point x. Not having independent evolution 
equations for the eigenvectors of S restricts our ability to interpret the geometry of the problem. 
The significance of (f> in this context is that when it takes the value(s) <j) = (ir) this means that 
u has aligned (anti-aligned) with an eigenvector of S, although which one we cannot say. The 
so-called 'stretching rate' a is an interesting quantity to consider (£ = u;/|u;|) 

Ul ■ (J U) • Su> 

a(x,t) = = (6) 

which has been related to the vorticity through an elegant Biot-Savart integral |3J| . For the 
three-dimensional Euler equations, the scalar vorticity u = \u)\ simply obeys 

Dt = ^ (7) 
The scalar a is obviously an estimate for an eigenvalue of the strain matrix S and lies within 
its spectrum. By defining the vector 

u> x Su> 

X(x,t) = =ZxS£ (8) 

the angle 0(x, t) between u and Su can naturally be introduced as (x • X = X 2 ) 

\u) X Su\ y , . 

and its dynamics related to those of a. and X- Gibbon and Heritage have recently discussed 
this angle but only in the context of the volume average over the whole flow. The stretching 
rate a takes account only of the symmetric part S of the deformation matrix Uij but the 
antisymmetric part clearly contributes to X- For a fluid element at a point x certain specific 
values of the angle <j> can be interpreted as the following: 
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1. When 4> = the vectors uj and Suj are parallel and only stretching occurs. 

2. When <f) = ir/2 the vectors u and Su are orthogonal and, since Su = Du/Dt, this means 
that in this case u rotates but does not stretch. This effect will nevertheless distend and 
misalign vortex lines. 

3. The case <p = ir represents anti- alignment. In this case a < 0, so vorticity collapses rapidly 
to small values. 

4. When < eft < ir/2 both stretching and rotation occur simultaneously at each point on a 
vortex line. 

In order to find the most sensitive relationship between a and x we take advantage in §2 of a 



result of Ohkitani ]4g, 47] for the incompressible three-dimensional Euler equations 

where P = {p,ij} is the Hessian matrix of the pressure. This result^] comes about from a more 
general well known property of the Euler equations, namely that if u and u are velocity and 
vorticity solutions of the Euler equations then for any arbitrary vector A 

dAi d DAi , s 

r-Wj- — (11) 

oxk oxj Dt 

from which we conclude that 

= ( B .VA)=«.V(^). (12) 
Choosing A = u with a direct use of the velocity form of the Euler equations on the right hand 




side of equation fll2[) immediately gives equation (10). While the pressure Hessian in (|i0| ) is 
a fully nonlinear term, (trP is related to Uj by Ap = —UijUjj), nevertheless the cancellation 
of two terms inherent in the derivation of equation (]Io|), each of which are ~ o^|Vu| 2 , directly 
removes all the terms not explicitly dependent on P. The approach of this paper therefore runs 
counter to the conventional one where the pressure is removed by projection but the advantage 
gained by the fortuitous cancellation of nonlinear terms comes at the price of having to deal with 
the P-matrix. Previous studies on depletion of nonlinearity have been made by Constantin and 



Fefferman [45] for the three-dimensional Navier-Stokes equations and Constantin, Fefferman and 



Majda [39] for the three-dimensional Euler equations who have considered the angle between 



the vorticity vectors u;(x) and <*»(y) at two points x and y in the flow. 

Note added in proof: In a private communication, P. Constantin has informed us that the result in equation 
(10) has been known privately in unpublished form since 1983. 
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The main results of the paper can be summarized as follows. It is shown in §2 that at points 

in the flow for which u ^ 0, the material derivatives of a and x f° r the Euler equations obey 

D ° 2 2 fia\ 

— = x - a -a p , (13) 
Dy 

-j£ = -2a X - X P , (14) 

where a p and Xp are related to the Hessian matrix P which is a non-local quantity in the 



flow jl7], 48, |4S| ]. For the Navier-Stokes equations (see §3) the addition of viscosity produces 
differential equations of the type 

= x 2 -a 2 + vka + 2vu\Vl\ 2 + \ (15) 

-^ = -2ax + ^A X + 2^|V£| 2 + //, (16) 
where the terms in |V£| 2 express the misalignment of vortex lines in the differential geometric 



sense [50, 51]. The terms A and fi contain information about the flow and are not generally 
constant. In §4, however, the exact solutions of the Navier-Stokes equations that represent 
Burgers vortex tubes and shear layers are discussed and it is discovered there that these corre- 
spond to 'Lagrangian fixed point' solutions of (|l!]) and ( |T6| ) in the sense that Da/Dt = and 
Dx/Dt = 0, A and a are positive real constants and = and X = 0. Hence for both these 
thin structures A and fi simplify greatly and the corresponding angle is exactly zero. Using 
this as motivation, it is assumed in §5 that if the flow is regular then all the variables must come 
to some equilibrium in a connected region of high vorticity. From this we want to see if the 
equilibrium values of a and x correspond to a small angular orientation <f)Q. In order to do this 
it is assumed that A and fi are constant and then it is shown that equations (|15|) and ( |i~6| ) have 
two fixed point solutions (cto,Xo) * n the Lagrangian sense, one of which is repelling (ao < 0) 
while the other (ao > 0) is attracting. The stability of the positive root for a shows that the 
system favours vortex stretching, in agreement with G. I. Taylor's conclusions [I]. For A > 0, the 
stable solution has a corresponding small attracting angle (j)Q which is insensitive to the relative 
ratio of A and making the natural orientation of the vectors close to true alignment. These 
results are consistent with near Burgers-like structures forming in the flow. It remains to be 
proved, however, that the Burgers solutions are either unique or belong a more general unique 
class of solutions that correspond to thin sets which are true Lagrangian fixed points with small 
values of eft. As we point out in §2 neither Hill's spherical vortex (an exact solution of the Euler 
equations) nor ABC flow belong to the class we are considering as neither are Lagrangian fixed 
points of their respective (a, x) equations. 
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2 Stretching equations for the three-dimensional Euler equa- 
tions 



The scalar a and the vector x were defined in equations (g) and (|8|) by forming the dot and 
cross products respectively of uj with Suj. We repeat their definitions: 

uj • Suj " " uj x Suj 

UJ ■ UJ • UJ 

but exclude points where u) = 0. The unit vector x(x, i), defined by 

* = -, (18) 

X 

is orthogonal to the unit vector for example, if £ points along a vortex tube then x nes i n the 
plane orthogonal to the tube. Let us also define the scalar quantity a p and the vector Xp 

, . us ■ Puj , . uj x Puj , , 

a p {x,t) = x(x,t) = . 19) 

UJ • U y U) ■ U) 

a p is an estimate for an eigenvalue of the Hessian of the pressure P = {p,ij} in the same way 

that a is for the strain matrix S. When u aligns with an eigenvector of P then a p becomes an 

exact eigenvalue of P whereas x P becomes zero in this case. 

In §1 the definitions @ and (||) were used to define the angle <f> between u and a = Su>, 

namely tan</> = x/ a i as a wa Y °f characterizing the alignment or misalignment between them. 

In terms of the vector x this angle obeys (<x = Suj) 

uj x a uj x Suj , . 

Xtan^= = —, 20 

ui ■ a uj ■ Suj 

with <j) varying in the range < <fi < 2ir. One consequence of equation (^0|) is that 

o , (er ■ uj) 2 uj 2 o? , . 

cos 2 ^)= , V .-. ; , 9 = — 21 
|cr|^ \uj\ l o A 

To obtain relations between these quantities we firstly note that 

D(uj ■ a 



Dt 

whereas in the cross product, the terms in a vanish leaving only 

D(uj x a 



<j — uj • Puj = uj ( of sec <j) — ot p j . (22) 



— uj x Puj. (23) 

Dt y ' 

Equation ( |23| ) shows that in addition to the cancellation of nonlinear terms in (10), more non- 
linearity has been lost through the cross product. The material derivatives of a and x can now 
easily be obtained through equations (||), (p|), ( p2[ ) and (23), 

Da 9 9 , 

— = x 2 ~ « 2 - op, (24) 
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-2a X - X P - (25) 



Dt 

Subject to solutions existing, the prominent feature of ( |24| ) and ( p5| ) is that they are independent 
explicitly of u and (f> while being driven by a p (x,t) and X p (x, t). Hence they are a set of non- 
autonomous ODEs, in the Lagrangian picture of fluid mechanics, operating in a 4-dimensional 
phase space (a,Xi) for i = 1,2,3. The three equations in the Xi can De reduced to one in the 
scalar x = \fX • X to give the pair of equations 

— = x 2 - « 2 - ap, (26) 
Dy 

-^ = -2a X -X P , (27) 

where 

X P = X ■ X P - (28) 

On the a-axis where (j) = or tt, equation (|2^) show that u x Pu = 0. Hence x P = an d so 
Dx/Dt = on this axis. Consequently equation ( p7| ) operates only in the the upper half-plane. 
The lower half-plane is simply a reflection of the upper in the a-axis. While a is the stretching 
rate for the scalar vorticity u>, equation (|27]) shows that 2a is the contraction rate for x with an 
additional pressure term. The two scalars u and x are also related by 

D 2 uj / o \ / % 

— = (x 2 -a p )u J . (29) 

The pair of equations for (a, x) given above can be merged into one by defining the complex 
stretching rate 

C = ol + i X (30) 

and (|2^) and (|?]) reduce to 

^| + C 2 + Cp = 0. (31) 

Cp is simply the complex combination £ p = a p + ix p . Because the C-plarie is only the upper half 
of the complex plane, <f> is valid only in the range < (p < tt. Two related results, which are 
consequences of equations (|26| ) and (p^) , are 

D(\C\ 2 ) 
Dt 

and 



-2a|C| 2 -23f?{cC p *}, (32) 



D{t ^ ] - -atan 3 0-fa-^) W-^- (33) 



Dt \ a J a 

Because of the negative cubic term, there is always the tendency in the a > quarter-plane for 
the angle to decrease unless the pressure terms force it to behave to the contrary. Equations (|2q) 
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and (|27D, however, are non-autonomous and, for Euler flows, there can be no question of a p and 
X p remaining constant as the flow develops. Nevertheless certain tendencies can be discussed: 

1. In the left hand quarter-plane equation where a < (|7|) shows that vorticity must be very 
small here. Initially if one starts in the region a < then equation ( |27| ) shows that x can 
undergo rapid growth with 4> — > tt/2. This is also reflected in equation ( |32|) and equation 
( |33| ) where the terms — 2a|£| 2 and —a tan 3 (f) force rapid growth in |£| and tan (/) respectively. 
If x increases in this way then ( p6[ ) shows that there is also the tendency for a to increase 
also. Consequently there is a natural tendency for orbits to pass from the left to the right 
hand half plane unless the pressure acts to prevent it. 

2. In the a > quarter-plane the opposite process occurs. Here, if a p < a 2 and x p > for a 
sufficiently long time then — > 0, which is the state of exact alignment. If a p > a 2 then 
how far 4> decreases depends upon specific values of a p and x.p- 

We conclude that the system has a natural tendency to reject negative values of a unless forced 
by the pressure Hessian to behave otherwise. This means that the system tends to prefer regions 
where stretching is positive. It is noteworthy that in rejecting the a < region an orbit must 
pass through the x _ax i s into the region of a > 0. This process can be reversed if the pressure 
behaves in a contrary manner but this behaviour can only be properly elucidated from data and 
this analysis is in progress. 

What of exact solutions of the Euler equations? As we shall see in §4, the Burgers vortex 



solution of the Euler equations is a Langrangian fixed point of pq ) and (27) with a = const, x = 
with a corresponding angle <j) = 0. Hill's spherical vortex, on the other hand, also has (ft = 
but this does not correspond to a Lagrangian fixed point as fluid packets do not have values of 
a for which Da/Dt = 0, based on the fact that u • Va 7^ 0. To see this consider a vorticity 
field u = (uj r ,ujg, uj z ) where co r = u> z = and ujq = Ar inside the sphere and zero outside it p2|. 
Then a = Az/5, x = 0, a p = A 2 Ur 2 - ^f-] /50 and x p = 0. The solution for ABC flow for 
which u = u with u\ = u\ = sinz + cosy, 112=^2 = sinx + cosz, 113 = L03 = s'my + cosx, 
does not have constant values for a or x- This means that neither cj) is constant nor are the 
material derivatives of a and x zero. 
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3 The stretching variables and the Navier-Stokes equations 



For the Navier-Stokes equations, the task is to see how the addition of viscosity changes the two 
equations for a and x in ®) an d (p5|). In vorticity form the Navier-Stokes equations are 



Dt 

and the equivalent of (0) for the scalar vorticity uj = \oj\ is 



- Sj k uj k + vAujj, (34) 



-j^j- = auj + v£ ■ Au>. (35) 



/,From the definition £ = lo/uj 



^ = S jk £ k - <• + -Acjj - -(£ ■ Aw)&. (36) 

Ut UJ UJ 

It is convenient to turn derivatives of cj into derivatives of To achieve this we use the fact 
that = 1 leading to CjA^j + |V£j| 2 = 0. From the relations u = uJj^j and ujj = u£j we have 

Auj = w|V£| 2 + £ • Aw, (37) 

and 

= A& + 2V (In u) ■ V& +6—- (38) 
uj u; 

Dividing ([}?]) by uj, multiplying by £j and then adding to (|38|) enables us to rewrite (|^) as 
where the vector /3 has components 



SjkZk - ocij + vA£j + i/|V€|%- + (39) 



^• = ^-(ln|a,| 2 )|^. (40) 



|V£| is the 'misalignment' of the vector £ in the Frenet-Serret frame (see references [50, and 
the calculation in Appendix 1). Equation^] ( |40| ) is the key to the calculations in the rest of this 
section because, apart from the logarithmic derivative in u, it contains only £ and its derivatives. 
The strain matrix evolves according to 

^ - -Sifc5 fcj - + - Pij + vASij, (41) 



and so, 



+ v { {ASitfj + SijAZj + | VtfSijtj + Sijfy } . (42) 



2 This expression is not valid at stagnation points where u> = neither is £ defined at these. 
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This, together with (|39|), can be combined to give, 

§i& S i&) = X 2 ~o? + 2v\Vl\ 2 a 



+ viZiSyPj + PiSytjy-tiPijtj. (43) 



X 2 - a 2 + v/\a + 2v\V£\ 2 a + ^A jj , (44) 



Out of the Laplacian terms in ( f43[ ) it is desirable to form Act 

Da 
~D~i 

where 

Aji = v (tiSijft + PiS^i) - 2uT fl - ^P -e,, (45) 
and where Tji is defined as 

with /3j is defined in (@Q). We have isolated in (44) as many specific terms in a as possible and 
in forming the Laplacian term Aa we have separated the highest derivatives of S from the rest 
which lie in Tj\. For x a similar result holds which is 

Dxi 



Dt 

In summary, if we define 



-2a Xi + v&Xi + 2HV£| 2 Xi + ^jfc^fci- ( 47 ) 



Hi = £ij k A kj , and A = ^A^-, (48) 



then equations ( f44| ) and fl47| ) become 

Da 



and 



_1 = x 2 -a 2 + i/Aa + 2^|v£| 2 a + A, (49) 



-^ = -2ax + z.Ax + 2^|V^| 2 x + /i. (50) 



As in the Euler equations in §2 this can be reduced to one equation in the scalar x 

^ = -2 X a + vA X + v (2\Vl\ 2 - |Vx| 2 ) X + A, (51) 

where 

M = X ■ A*- (52) 
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4 Burgers vortex and shear layer solutions 



We now postpone our discussion of the (a, x) equations for the Navier-Stokes case until §5 
and return to the Burgers vortex and stretched shear layer solutions. It is pertinent to ask 
whether equivalent Langrangian fixed point solutions exist for both the Euler and Navier-Stokes 
equations? Despite the caveats made in §1, the simplicity of the Burgers solutions makes them 
candidates for this. In the following two subsections the basic formulae for the axisymmetric 
Burgers vortex^] (see p5| , 13, ^4], and the Burgers shear layer (see |l3|, |24|, |35|) are worked 



out. Then in §4.3 we see how these are applied to the (a, x) equations ([49]) and (^) for the 
Navier-Stokes equations and find that they do indeed correspond to fixed point solutions. In 
the paper by Majda a series of simple examples is given which includes pure rotation, a 
swirling drain and the two Burgers solutions (see §1C of reference [|l3|]). The two expressions for 
the vorticity ^>(x, t) for the latter pair of examples are based on an ^-dimensional heat kernel 
but for the sake of simplicity in the following two subsections we use the time asymptotic limit 
of this. 

4.1 Burgers vortices 

Consider a strain field u = (— — ^■, r yz) T with a superimposed two dimensional velocity field 
(— y/(r), xf(r), 0) T where the function f(r) will be fixed later on. The variable r is given by 
r 2 = x 2 + y 2 . Then the full velocity field is 

u = f-y - Vf(r), - y + xf(r), -yz\ (53) 
and it is easily seen that the vorticity is 

u = (0,0,u, 3 ) T , (54) 
where uj^ = 2/ + rf. Moreover the strain matrix S is given by 



/ 7 xyf (x 2 -y 2 )f n \ 

2 r 2r u 

^ 2 -v 2 )f _z + mL n 

2r 2 ' r u 



(55) 



V 7 J 

Immediately we see that e 3 = (0, 0, 1) T and A3 = 7 and the other two eigenvalues can easily be 
computed 

Ai, 2 = (56) 



3 Note added in proof: See also the paper by Andreotti B. (1997) "Studying Burgers' models to investigate the 
physical meaning of the alignments observed in turbulence," Phys. Fluids A 9, 735-42. 
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with their corresponding eigenvectors lying in the horizontal plane. Now for the Euler equations 
one can simply take f = -K Jq sus(s) ds and so u = (0, 0, u^(r)) T . For instance, there is a 
solution of the three-dimensional Euler equations on compact support whose amplitude increases 



exponentially while its support decreases exponentially [19]. 

For the Navier-Stokes equations, however, there exists a form of f(r) in which dissipation 



and stretching balance which, in the limit t — > oo, is given by [13] 

1 — exp (— ar 2 ) 



-Uq. 



(57) 



where a = ^y/Au. This is a profile which has its maximum when r = 0. Note that when r = 
then /(0) = a, u}^{&) = 2a whereas [r/'(r)] r=0 = 0. Near r = 0, rf ~ —2a 2 r 2 < and so we 
have two negative eigenvalues and one positive for this form of solution. Also the total strain at 
r = is 

37 2 
2 



(58) 



Hence if v <C 1 then a»7 and so the vorticity ^(O) = 2a is much larger than the total strain. 
An x — y plane contrast cross-section of the vorticity field for a Burgers vortex is produced in 
figure la. 

4.2 The Burgers shear layer solution 

Consider a jet u = (0, — 7y, Jz) T which compresses in the y-direction but expands in the z- 
direction. Now impose a velocity field v(y) on the x-direction so that 

u= (v(y),-7y,jz) T . (59) 

u = (0,0,u; 3 ) T , (60) 
where uj 3 = —v'(y). For the Navier-Stokes equations, in the limit t — ► oo, if v(y) is taken to be 

ryV2a 



1 rw/2a 

-= I v(0)exp(-s )ds 

V 27T J — oo 



v(v) 

with a = j/4u then 

CJ3 = \/ — exp(— 2ay 2 ) / wo(s)ds. 

V 7T 7 

This the Burgers shear layer for which the strain matrix S is given by 

' y N 

v 7 J 



(61) 



(62) 



(63) 
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As before we see that e3 = (0, 0, 1) T and A3 = 7 and the other two eigenvalues are 

^ ^'^^ (64) 

with their corresponding eigenvectors lying in the horizontal plane. In this case the total strain 
is 

E4 = 2 7 2 + ^- (65) 

When v -C 1 then a S> 7 and therefore |t>'| 2 3> 7 2 . Hence when the vorticity is high then so 
is the total strain, in contrast with the vortex tube in the previous subsection. An x — y plane 
contrast cross-section of the vorticity field for a Burgers stretched shear layer is produced in 
figure lb, 

4.3 The (a, x) equations for both tube and shear layer 

Taking the third component of the velocity version of the Navier-Stokes equations it is easy to 
show in both the above cases 

I = -7* (66) 

and hence 

a p = - 7 2 , X P = 0. (67) 

In addition, using the solution for u and the form of S, in both cases 

a = 7, X = 0, (68) 

with corresponding angle 

<t> = 0, (69) 

reflecting the fact that there is exact alignment between u> and e3 with this fixed point sitting 
on the a-axis. Clearly, in terms of the angle we can make no distinction between the tube and 
the shear layer, nor is there any distinction between the cases of alignment or anti-alignment of 
uj and an eigenvector of S. Moreover, the unit vectors £ and x are given by £ = k and x = so 
V£ = 0. This means that when fx = and A = 7 2 , the relevant (a,x) equations given in ( |4^ ) 
and (|5^) are satisfied by equations (^) and (|68|) for which 

Da Dx , x 

The problem of stability is considered in the next section; there we will show that this fixed 
point is indeed stable provided A and fi remain constant. Given the fact that both Burgers cases 
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are essentially 2D velocity fields superimposed on three-dimensional strain fields, any example 
of this type will produce a strain matrix S with a similar block diagonal form as in ( |55| ) and ( |63| ) 
with a corresponding eigenvector e3 which is parallel to the vorticity vector u. All examples 
of this type therefore have = 0. Moffatt, Kida and Ohkitani |gq] have used asymptotics on 
the Burgers vortex where they allow a certain small asymmetry from the axisymmetric solution 
(E7j) but they point out that only the exact symmetric solution is known. 

What is not proved, however, is that if the system is attracted to a stable Lagrangian fixed 
point solution with an orientation near eft = then this must automatically correspond to either 
of the Burgers-like solutions sketched above; there may be other unknown structures which fall 
into this category (see the comments at the end of §2 on Hill's spherical vortex). The lack 
of certainty in the topology reflects the problem of not having enough dynamic angles in the 
system. 

5 A mechanism for the formation of thin structures in Navier- 
Stokes flows 

So far we have made no assumptions. Let us now discuss some approximations and put forward 
the following theoretical picture which is consistent with the formation of thin structures in 
isotropic Navier-Stokes turbulence. Over all calculations of this type lies the heavy shadow of 
the technical question of regularity, which is still an open problem. The fact that alignment in 
mature turbulent flows is being discussed, however, means that a sufficient degree of regularity 
is being imputed to the solutions anyway. Prom now on we assume that the flow is regular and 
that all the necessary quantities are bounded. ^From §4 we know that when the flow assumes a 
Burgers structure, A and /x take the values 

A = 7 2 , ti = 0. (71) 

In the following, however, we do not assume that A and fi take the values given in ( |7l|) as this is 
tantamount to assuming the answer. Nevertheless, the highly simplified form that A and fi take 
when the flow assumes a Burgers state motivates us to assume that the variables a, /x, A and 
\x have reached a simultaneous equilibrium in some connected region of vorticity whose growth 
has been controlled by dissipation. First we want to see if the 'fixed point' values taken by a 
and x a t this equilibrium correspond to a small angle (fio and secondly whether this fixed point 
is stable. To investigate this stability question we begin by assuming that A and /x are constant 
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in the four equations in fl49|) and d5_ 



Da 
~Dt 



X 



a 2 + uAa + A, 



— = -2ax + z^A% + n 



(72) 



(73) 



where, for reasons of simplicity, we have written these without the quantity |V£| 2 . These terms 
are the least important and are dealt with in the Appendix later. Fixed points in (ao, Xo) occur 
at 

" (74) 



2ai 



A + V V + A 2 , 



Xifi 



JM_ 

2«o 



where /x = \fjt\. Hence there are two fixed points in 4-space, corresponding to the two roots of 
ao in (|74j). Without the Laplacian term it is easy to show that the eigenvalue stability problem 
is 



(A + 2a ) 2 [(A + 2aor + 4xj) 

thereby giving the four roots 

A = — 2ao (twice), 



0. 



A = -2 (a ± ixo) 



(75) 



(76) 



These roots correspond to an unstable fixed point for the negative root for ao and a stable one 
for the positive root. In the latter case, there is exponential contraction in two of the directions 
in the 4-space with a stable spiral in the other two directions. When the Laplacian is included 
we look at the stability of linearized solutions of the type exp (ik ■ x + At) around (ao, Xo)- 
only difference this makes to fl76|) is that 2ao — ► 2ao + vh 2 . We conclude that provided ao > 
the equilibrium solution is stable to the disturbance of all wavenumbers. For the half-plane the 
two fixed points are 



A + vV + A 2 



ql/2 



V2 



Xo 



A + vV + A 2 



■lV2 



(77) 



The stable spiral is shown in figure 2 which is projection onto the variables xi an d a. Defining 
m to be 

f 1 



m 



|A|' 



(78) 



the angle (f>o corresponding to the stable fixed point is (the + sign is for A > and the — sign 

for A < 0) 

. _ \/m 2 + 1 =F 1 



tan < 



771 



which can be simplified to 



tan 2(/>q = ±7ii. 



(79) 
(80) 
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Now the exact Burgers solutions both correspond to m = and, because a = 7 > 0, they cor- 
respond to an attracting fixed point of the system. We note however that when it takes nonzero 
values, m is dependent only on the ratio of \i and |A|. Although we have no hard information on 
the magnitude of m, in fact the angle </>g is relatively insensitive to this magnitude. For A > 0, 
if m ~ then </>q~ ~ but if m ~ 1 then </>q~ ~ 7r/8. Even if m = 00 then (pQ = 7r/4. Hence 
at worst, lies in a 45° cone. Therefore, even when m > we are still close to alignment. 
More generally, /x and A derive respectively from forming the vector and scalar products of the 
same set of functions so, in an isolated region of vorticity f2, we might expect that on the spatial 
average over a random set of points in Q, m would take the value m ~ 1. It is also possible 
that because /x is formed from a cross product while A is formed from a dot product then at 
the natural angle of alignment the \x term would be the weaker of the two making m < 1. We 
therefore conclude that over the spatial average within a small intense region for A > 

0+ « tt/8 (81) 



but that 4>q may indeed be somewhat smaller than this. Experiments [ESJ, 30, BlL p2j and 



simulations (LJ, 16, 18, 2C, |28[ generally measure the cosine of the angle <f>, often observing 
that a bunching around cos (j) ~ 1 in their PDF's is a demonstration of alignment. For <^>q" ~ tt/8 
we have 

cos #f ps 0.92. (82) 

This value of m therefore produces fairly close alignment. This is consistent with the region 
finding an equilibrium shape near to a thin Burgers-like structure (for which the exact value of 
m is zero) which corresponds to A > 0. When A < 0, however, 

<P = vr/2 - <f>+ (83) 

and the equivalent value of <pQ is 3ir/8. Hence vortex lines are badly misaligned in this case. 



6 Conclusion 

The thin vortical structures which are observed in numerical simulations of turbulence have 
been one of the most intriguing visual manifestations of the complexity for which turbulence, 
for better or for worse, has become a byword in the last two decades. Our results are consistent 
with near Burgers-like structures (see figure 1) forming out of those connected parts of the flow 
where A and /x equilibrate and where A > 0. The fact that the stable equilibrium occurs in a 
region of qq > means that there is a preference for vortex stretching [JiJ. To say that these 
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regions consist of precisely tubes and/or sheets is to oversimply the matter; more generally the 
individual structures probably have a fractal dimension which lies somewhere between one and 
two corresponding to tubes and between two and three corresponding to sheets |2], [3|, [|, p3|] . 
Nevertheless, tubes and sheets as extreme limits of the topology suffice as an approximate 
description despite the fact that their local interactions and constant metamorphosis produces 
levels of geometric complexity which are beyond our understanding at the moment. This raises 
several questions regarding the geometric consequences of vorticity alignment: 

1. Can all solutions for which eft ~ 0, which are also Lagrangian fixed points of the (a,x) 
equations, be categorized as 'thin' in the sense that they are quasi ID or 2D? While we 
have characterized these structures by the angle </>o(x, i) which lies between to and Su, are 
there other dynamic angles which characterize the topology more specifically? 

2. Are the Burgers solutions or their generalizations unique among the set of thin solutions? 

3. What is the behaviour of m = /u/|A| defined in ( |78|) and what determines the sign of A? The 
assumption that A and n are roughly constant in some regions is based on the assumption 
that (consistent with the fact that they behave this way for the two Burgers solutions) a 
balance occurs at the deepest scales between the pressure and the viscous terms expressed 
through P, S, VS 1 and Vw. Such an assumption would be hard to prove directly while no 
regularity proof exists. 

4. How are a and x affected when there are local interactions; for instance, when two sheets 



interact and wrap up to become tubes [26]? It is possible that large values of A correspond 
more to sheets than tubes since the former have higher strain. The experiments of Andreotti 
et al 1 36] creating Burgers vortices suggests that the tubes formed from the roll- up of sheets 
have a long lifetime compared to the roll-up time. This suggests that two time scales may 
be involved; the first being the time it takes for the orbit in (a, x) space to stabilize and 
the second the time over which A and fi remain constant before themselves decaying due 
to dissipative effects. 

5. What role does the Hessian matrix P play in these processes? Is there a consistent pattern 
in the sign of a p ? 

As far as question 1 is concerned, it is clear that one angle is by no means sufficient to adequately 
describe the topology. There are two other apparent natural angles in the system. The first of 
these is based on the matrix Uij instead of which we discuss only in the Euler case. Define 

a H = = Yh ( 84 ) 

U ■ UJ UJ ■ UJ 
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and then use the general result ( |12[) to show that 

— — — = u> ■ Va. (85) 
Dt y J 

~~ D shows that points moving with the flow in regions where a = const (or even just spatially 



(86) 
(87) 



independent) also have a^ 3 ' as a constant of the motion. Now define 



UJ X {uJ ■ V)u? 



then 



— = w • V X + 2a^x - 2a X (w) + * ^. 

Dt uj ■ uj 



When (ft = then er = = auj and X = and so 



0. 



In consequence the angle ^ between a> and (a; • V)a; for points travelling with the flow is given 
by 

const , . 

tan-(/>=^^. (89) 

tp therefore decreases in regions where aM is stretching. The second angle concerns the Hessian 
matrix of the pressure P 

\uj x Poj\ , . 

tan</» P = J —i-, 90 

CJ • Puj 

but we have no separate knowledge of the evolution of P. In fact P and S are less independent 
than one would think. A consequence of equation (|2^) is that when uj aligns with an eigenvector 
of S then uj also aligns simultaneously with an eigenvector of P. 

The lack of knowledge of m raised in question 3 above is balanced in part by the insensitivity 
of 0o, the attracting angle, to the value of m through the relation (A > 0) 

tan20 o f = m. (91) 

There is no reason, a priori, why exact alignment at m = should be attracting. Nevertheless, 
values of m > produce attracting orientations which are still close enough to alignment to 
suggest that we are near a thin structure which we conjecture is perhaps twisted in some way. 
As we showed in the last section, different regions may take different spatial values of m but 
their associated angles may still be close enough to look the same in observations or simulations. 
m obviously contains geometrical information about the local and relative orientation of vortex 



lines in the sense of the Frenet-Serret equations [ 50 , 51 1 but these need information on Va; which 



we do not have. It is possible that more can be said about m using a scaling argument. 
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What the formulation of this paper does not do is differentiate between tubes and sheets; the 
alignment process is the same for both and we have no other information which distinguishes 
them. In fact, one very important set of processes, not accounted for by the dynamics of 4> alone, 
are the interaction processes between one sheet and another. The Kelvin-Helmholz instability, 
which is an Euler phenomenon, is a well known mechanism through which it is thought that 



two dimensional sheet-like objects wrap up to become tubes [38, 27, 20, 54|. Such local 
interaction processes are extremely subtle and if they are to be accounted for dynamically then 
certainly more than one angle would be needed. 

The thin structures observed in simulations take up no more than a few percent of the total 
flow volume indicating that vorticity is not distributed evenly. This may be one of the reasons 
behind the failure to prove regularity. The methods used have been based upon attempts to 
control the global enstrophy, Jq \oj\ 2 dV, which have foundered on the problem of the dissipation 
being too weak to control the estimated nonlinear terms (see references in |55|]). These latter 
terms have their origin in the vortex stretching term, J ojiSijLOj dV, and have to be estimated 
by standard Sobolev inequalities. The failure to control the global enstrophy for more than a 
finite time by this method probably has its origins in the use of volume integration in which 
L 2 -norms average over, and perhaps miss, the spatial structures where the nonlinearity has been 
strongly depleted by alignment. Constantin's Biot-Savart integral formula for a in terms of a 
triad of vectors related to £ illustrates how nonlinearity is depleted inside the volume integral 
when significant alignment takes place [f40| . 

The general picture that emerges from this analysis is that the variables (a, x) seem to form 
a natural pair of variables in which to express the dynamics of vortex formations. In both 
the Euler and Navier-Stokes cases negative values of a are repelling and, without interference 
from the pressure, <j) would limit to zero. Any initial state corresponding to a negative value 
of a would require the system to go through the <j) = ir/2 stage before becoming a Burgers- 
like shear layer or tube if indeed it is driven that far. The difference between the two is that 
for the Euler equations there is no stable equilibrium point provided by the viscosity and the 



whole process could reverse, although Ohkitani and Kishiba [47] report that they observe an 
alignment between u and the third eigenvector of P. A systematic study of simulation data 
for the three-dimensional Euler equations is underway to understand the properties of the more 
general behaviour of the Hessian matrix P. 

One curious conclusion from (^6|) in the Euler equation case, for a Lagrangian particle element 
initially at X = x(0), is that if a(X, 0) < and </>(X, t) < 7r/4 with a p > then a — > —oo in a 
finite time. In contrast, if a(X, 0) > and </>(X,i) > 7r/4 with a p < then a — > oo in a finite 
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time. This is consistent with the results of Ng and Battercharjee who have studied the Euler 
equations under high symmetry conditions which makes the problem quasi one-dimensional. 

Past simulations have shown that initially cj aligns with the largest eigenvector of S but 
as the turbulence becomes more mature u aligns with the second eigenvector []l4| , [l6| , |26| as 
a global average. Tsinober, Shtilman and Vaisburd |?2| have recently reported, however, that 
locally they see significant alignments between u and both the first and the second eigenvectors 
of S. They also report that the vorticity does not have to be too intense and they produce 
evidence to suggest that the background field is not Gaussian. 
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A The effect of the misalignment terms 

For the | V£| 2 terms left out of the discussion in §5, it would be preferable if we had five evolution 
equations in a, x% an d V£ and not just four but we have been unable to discover any delicate 
cancellations in the material derivative of the last variable. Instead, we include this as a constant 
term in the equilibrium point analysis above. For a vortex line, V£ is called the misalignment 
and plays an interesting role. Constantin, Procaccia and Segcl |5(J note that 

|V£| 2 = \(i ■ V)|| 2 + |(n • V)£| 2 + |(b • V)£| 2 (92) 



where the first term on the RHS of (92) is the square of the curvature and the second and third 
are the squares of the the lack of parallelity between a vortex line and neighbouring lines in 
the normal (n) and binormal (b) directions respectively. Galanti, Procaccia and Segel [|5l] have 
produced theoretical and numerical arguments which indicate that the stronger the vorticity 
and the greater the curvature of a vortex line the stronger the stretching and therefore the more 
liable it is to straighten. To show how this is the case we define 

a = zv|vt| 2 a = a-a A = A + a 2 . (93) 

The (a, x) equations become 

^ = X 2 "« 2 + A, (94) 
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Dy 

j± = -2a X + fi. (95) 



The stability problem is unchanged and fixed points come from 

1/2 



A + (\ 2 + fi 2 ) 

4 = [ o J ■ ( 96 ) 



^From (96), Oq > A > a when A > and so ao > 2a for positive values of ckq. Moreover, when 



considering the angle of orientation of the equilibrium point 0o we see that 

\jA = M 

|A| |A + a : 



tan 20o = H = 1^1 < m . (97) 



Hence 

tan 20o < tan 2</>o . (98) 

Equation ( p8| ) shows that nonzero values of |V£| 2 act to make the angle smaller than when it is 
excluded. Of course the angle and |V£| 2 are different ways of expressing the same effect. We 
actually need another differential equation to form an accurate picture and this, in turn, needs 
information on VS and Vu. The above calculation shows that for the (a,x) equations alone 



large values of |V£| force stronger alignment in 0, thereby confirming the effect seen in [51| 
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FIGURE CAPTIONS 



Figure la: x — y cross section of a Burgers vortex with the z-co-ordinate pointing out of the 
paper in which direction points the vorticity vector cj = (0,0,u>^) T . In the black and white 
contrast, the whiter the area the higher the vorticity. 

Figure lb: x — y cross section of a Burgers stretched shear layer. As in Figure la, the z-co 
ordinate points out of the paper in which direction points the vorticity vector uj = (0, 0,cu3) T . 
In the black and white contrast, the whiter the area the higher the vorticity. 

Figure 2: Orbits in the a — x phase space projected onto the a — xi plane. The spiral structure 
at the two fixed points is just evident. 
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